Galactic Halo substructure in the Sloan Digital Sky Survey: 
the ancient tidal stream from the Sagittarius dwarf galaxy 
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ABSTRACT 

Two studies have recently reported the discovery of pronounced Halo substructure in the 
Sloan Digital Sky Survey (SDSS) commissioning data. Here we show that this Halo substructure is 
almost in its entirety due to the expected tidal stream torn off the Sagittarius dwarf galaxy during 
the course of its many close encounters with the Milky Way. This interpretation makes strong 
predictions on the kinematics and distances of these stream stars. Comparison of the structure in 
old horizontal branch stars, detected by the SDSS team, with the carbon star structure discovered 
in our own survey, indicates that this halo stream is of comparable age to the Milky Way. It would 
appear that the Milky Way and the Sagittarius dwarf galaxy have been a strongly interacting 
system for most of their existence. Once complete, the SDSS will provide a unique dataset with 
which to constrain the dynamical evolution of the Sagittarius dwarf galaxy, it will also strongly 
constrain the mass distribution of the outer Milky Way. 



Subject headings: halo - 
kinematics and dynamics 

Introduction 



Galaxy: structure — Galaxy: halo — Galaxy: Local Group — galaxies: 
- galaxies: individual (Sagittarius dwarf galaxy) — galaxies 



In the standard hierarchical picture of galaxy 
formation, galaxies like the Milky Way grew by 
repeated mergers with other galaxies, and in par- 
ticular with many smaller galaxy "fragments" . 
In the outer regions of galactic halos, disruption 
timescales can be long, and we can hope to observe 
galaxy remnants as well as to catch galaxies in the 
process of being accreted. That this accretion is 
still an on-going process is clear, as evidenced by 
the many studies of tidal features in the dwarf 
satellite galaxies and globular clusters. 

However, what is of prime interest is to compare 
the observed structures to the results of galaxy 
formation theories. These now make strong pre- 
dictions about the structure and substructure of 
galaxy halos. Indeed the standard cold dark mat- 



ter models appear to be inconsistent with ob- 
servational constraints, as they over-predict halo 
dumpiness (see, e.g., Moore et. al 1999). 

In an attempt to constrain such theories, sev- 
eral groups have undertaken large scale surveys to 
uncover substructure in the halo. These include, 
among others, the "spaghetti" survey (Morrison 
ct al. 2000), the APM carbon star survey (Totten 
& Irwin 1998; Totten, Irwin & Whitelock 2000; 
Ibata et al. 2000), and the SDSS survey (Yanny et 
al. 2000; Ivezic et al. 2000). 

The APM carbon star survey, an almost all 
high-latitude sky survey, revealed the existence 
of a gigantic band of these intermediate age (~ 
6Gyr) stars around the sky (Ibata et al. 2000; 
hereafter paper 1). 4 Approximately half of the 
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4 Extant optical and NIR photometry and detailed optical 
spectroscopy of the APM colour-selected cool Halo carbon 
stars show that these stars are identical in nature to those 
found in AGB populations in other Local Group galaxies, 
particuarly the Magellanic Clouds and many of the dwarf 
spheroidal Galactic satellites (Totten & Irwin 1998; Totten, 
Irwin & Whitelock 2000). Furthermore the central part 
of Sgr has a known AGB carbon star population (White- 
lock, Irwin & Catchpole 1996) that also has NIR and opti- 
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"Halo" carbon stars detected in this survey are 
located in this ten degree wide band, which is 
centered on the expected orbit of the Sagittarius 
dwarf galaxy in a spherical Halo potential (pa- 
per 1). Sagittarius, the closest dwarf satellite 
galaxy to the Milky Way (Ibata, Gilmore & Ir- 
win 1994), is clearly in the process of being tidally 
torn apart. As shown in paper 1, those carbon 
stars trace a gigantic tidal stream of debris dis- 
rupted off the Sagittarius dwarf galaxy which has 
now wrapped itself around the Halo. Interest- 
ingly, this stream is observed to follow a great cir- 
cle on the sky, which means that it is probably 
a coplanar rosette-shaped structure that we ob- 
serve from within (since it should trace relatively 
closely the orbit of the center of mass of the dwarf 
galaxy). Most importantly, however, the fact that 
the stream is confined to a plane, implies that it 
does not suffer significant precession, which is only 
possible if the Galactic halo potential is not sig- 
nificantly oblate. 

Two SDSS teams have recently presented re- 
sults from their first year's commissioning data 
(Yanny et al. 2000; Ivezic et al. 2000). Due to ini- 
tial technical constraints, the drift-scanning cam- 
era took data only on the celestial equator dur- 
ing the commissioning period. Their data con- 
sists of two wedges: a NGC sample 87° long and a 
SGC sample 60° long, with each covering a decli- 
nation range of only —1.25° S ^ 1.25°. Despite 
the fact that their data provide essentially a two 
dimensional slice through the Galaxy, they find 
astonishing sub-structure in the distribution of A- 
stars. Sharply-defined concentric rings of stars are 
detected. The number of RR-lyrae variables also 
cuts off sharply at 50kpc. The nature of these 
structures is important, as they are the dominant 
structure seen in old Halo stars, comprising ap- 
proximately half of the total number of Halo stars 
seen in the survey. In this contribution, we show 
that these A-star Halo sub-structures are morpho- 
logically very similar to the expected structure of 
the tidal debris of the Sagittarius dwarf galaxy, 

cal properties indistinguishable from the APM Halo sam- 
ple. Normal AGB carbon stars are carbon-enriched during 
dredgc-ups driven by helium burning thermal pulses during 
the post main sequence phase of 1-2 solar mass stars. As 
such they are thought to be > 1 — 2 Gyrs old but less than 
10 Gyrs old (see e.g. the reviews by Da Costa 1997 and 
also Azzopardi & Lequeux 1992). 



are at similar distances and are highly likely to be 
part of the Sagittarius stream. 

At first sight it may seem surprising that such 
a limited survey can detect Halo sub-structure. 
However, since most long-lived Galactic satellite 
systems appear to be on close to polar orbits 
the chances of a stellar debris stream intersecting 
the current SDSS zero-declination survey strips is 
«50%. The Sagittarius stellar stream is currently 
the only known large relatively intact debris sys- 
tem and is within 13° of a polar orbit. 

2. The expected debris of the Sagittarius 
dwarf galaxy 

In paper 1, we undertook a large suite of nu- 
merical simulations of the interaction between the 
Sagittarius dwarf galaxy and the Milky Way, in or- 
der to determine how the APM carbon star sample 
constrains the shape of the Galactic halo. Fig- 
ure 1 shows a typical high resolution simulation 
with 10 5 particles, at the end of the integration, 
at T = 12Gyr (details of these simulations are 
given in paper 1). For this particular simulation, 
the initial structure of the dwarf galaxy is a fluffy 
King model with mass 1O 8 M , half-mass radius 
of r 1 / 2 = 0.5 kpc and initial concentration c = 0.5. 
The Milky Way model (termed model HI in pa- 
per 1) has a bulge, a thin disk, a thick disk and an 
interstellar medium, with parameters identical to 
the Dehnen & Binney (1998) model 2; the Galac- 
tic halo has a flat rotation curve, with a constant 
density core of radius r = 3 kpc, and an oblate 
spheroidal density structure, with density flatten- 
ing of q m = 0.9. The total mass of the Galaxy 
model is such that the circular velocity at 50 kpc 
is v c — 220 km s -1 . 

After accounting for a ~ 50% "background" 
from other accretion events, such simulations can 
give a reasonable representation of the observed 
carbon star distribution, both in spatial posi- 
tion and in kinematics (cf. paper 1). The Aitoff 
projection of this simulation onto the sky (Fig- 
ure 1), shows that the majority of the particles 
in this model are clustered close to a very obvi- 
ous great circle stream. However, the potential is 
not quite spherical, so orbits suffer some preces- 
sion. The material first stripped from the dwarf 
galaxy has had longest to drift away from its pro- 
genitor, and so also experiences the largest dif- 
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ference in precession. The particles not on the 
great stream are those that were stripped from the 
dwarf galaxy at earlier times; indeed, the preces- 
sion away from the great circle stream effectively 
provides a chronometer of the time since particles 
were disrupted from their parent galaxy. 

To compare the simulation of Figure 1 to the 
SDSS commissioning run data, we take a 2.5° 
equatorial cut though the simulation. As dis- 
cussed by Yanny et al. (2000) and Ivezic et al. 
(2000), there may be some mixture of blue star 
populations in their dataset, however, for simplic- 
ity, we assume that they are all blue horizontal 
branch stars with absolute magnitude M g * =1.0. 
The resulting distribution of these equatorial stars 
as a function of right ascension and apparent g* 
magnitude is displayed in Figure 2. The resem- 
blance of this distribution with that found by 
Yanny et al. (2000) (their Figures 3 and 11) is 
striking. This simulation (which preceded the 
SDSS result) predicts the azimuthal rings seen by 
the Sloan Survey. Figure 3 is identical to Figure 2, 
but the simulation it is derived from was under- 
taken in a more flattened halo, with q m = 0.5. 
The sharp ring features seen in Figure 2 and in 
the SDSS dataset are smeared-out in Figure 3. A 
spherical Halo (q m — 1.0) gives rise to a slower pre- 
cession rate (there is still some precession due to 
the modelled disk component), so a smaller frac- 
tion of particles wander off the great circle stream. 
Unfortunately, we cannot be quantitative in our 
comparisons to the models without a published 
list of SDSS A-stars. 

The color-coding of Figures 2 and 3 indicates 
the time elapsed since the particles left the tidal 
radius of their progenitor. In Figure 2, the tight 
clumps seen near 2 hours and 14 hours are due 
to recently disrupted material, which can be com- 
posed of young as well as old stars. This is con- 
sistent with the detection of a clump of carbon 
stars in these regions (paper 1). The particles that 
were disrupted long ago, displayed here in yellow 
and red, have had many Gyrs to precess away from 
the great circle stream. The only remaining visible 
tracers of this material will be old stars, such as 
RRLyrae variables. This model suggests that the 
ancient disruption and star-formation history of 
the Sagittarius dwarf may be traced and analysed 
with the population of RRLyrae variables first de- 
tected by SDSS. In the main body of the Sagit- 



tarius dwarf galaxy, blue horizontal branch stars 
make up a much higher fraction of the total stellar 
population in its globular clusters than in its field 
stars (Layden & Sarajedini 2000). The fact that 
the older parts of its tidal stream display appar- 
ently high concentrations of RRLyrae stars sug- 
gests a common origin with the Sagittarius dwarf's 
globular clusters (though this may just reflect con- 
tinual star- formation in the dwarf galaxy). 

Our simulations make strong predictions. First, 
the radial velocity of the A-stars should vary 
smoothly as a function of right ascension, as dis- 
played on the top panel of Figure 1. The radial 
velocity dispersion of these stars in a small range in 
right ascension should be small also, ~ 20kms _1 , 
except in regions where multiple streams cross 
each other, at a ~ 2 hours and a ~ 14 hours. 
Furthermore, the majority of the A-stars should 
have similar radial velocities to those of the carbon 
stars in the 2 hour and 14 hour regions of paper 1. 
On the equatorial plane of the sky (where the ex- 
tant SDSS data lie) at a ~ 2 hours, the velocity 
distribution should be bimodal, with two clumps 
at Heliocentric radial velocity of v ~ +200kms _1 
and v ~ — 150 km s -1 . Further along the equato- 
rial plane at a ~ 14 hours the mean Heliocentric 
velocity of the stream stars should be ~ 50 km s~ 1 , 
with significant dispersion of a ~ GOkms" 1 . Sec- 
ondly, the distances of the A-stars, and RR-lyraes, 
should also closely follow the range of distances 
found from the C-star studies. In particular, an 
upper cutoff of «50 kpc in the RR-lyrae distances 
is exactly what is seen in the C-star distribution. 

3. Conclusions 

We have shown that numerical simulations of 
the dynamical evolution of the Sagittarius dwarf 
galaxy, previously reported in paper 1, predict the 
presence of sharp shell-like features in the Galactic 
halo. This shell-like structure is due to the effect 
of orbital precession acting on stars which were 
disrupted from the Sagittarius dwarf galaxy many 
billion years ago. Making a cut through our simu- 
lated dataset along the equatorial plane of the sky, 
we find a distribution very similar to that found 
in A-stars in the SDSS data (Yanny et al. 2000; 
Ivezic et al. 2000). 

The A-stars furthest away from the dominant 
stream detected in paper 1, belong to those popu- 
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lations that were first disrupted from the dwarf 
galaxy. By comparing the relative densities of 
young to old stars in these tidal structures, it may 
be possible to reconstruct the dynamical evolution 
of the Sagittarius dwarf galaxy stream in great de- 
tail (Johnston et al. 1999). 

The SDSS data that will soon become available 
along other scan-lines on the sky offer the very 
useful opportunity to distinguish between the dif- 
ferent dynamical models of the Sagittarius dwarf 
galaxy. Such data will also allow us to distinguish 
competing models for the mass distribution of the 
dark halo. For instance, a comparison of the par- 
ticle distributions of Figures 2 and 3, shows that 
the enhanced stream precession in the flattened 
halo of Figure 3 gives rise to a substantial den- 
sity between 270° < a < 0°, which is not seen in 
Figure 2. Data in that right ascension range will 
provide a means to further constrain the shape of 
the Galactic halo. 

With the almost all high-latitude sky cover- 
age of the APM carbon star survey, only two 
Halo streams were found (the other probably being 
due to the Large Magellanic Cloud). Given that 
SDSS has also intersected the tidal stream from 
the Sagittarius dwarf galaxy, this stream appears 
to be by far the dominant contributor to Halo 
substructure at Galactocentric distances less than 
<~ 50kpc. It seems unlikely, therefore, that any 
more accreted systems of such magnitude will be 
found by SDSS or other surveys. However, there is 
hope to discover the much smaller accreted galaxy 
fragments, if they exist. 

The simulation result of Figure 1, has a cau- 
tionary tale to tell about using starcounts to de- 
termine the shape of the "Halo" . That a substan- 
tial fraction of the outer Galactic halo stars ap- 
pear to be clumped on an inclined band on the 
sky clearly does not mean that the Halo is a thin 
inclined flattened band! Such structures just tell a 
complicated story about the process of stellar ac- 
cretion into the Halo. Since stars only contribute 
a miniscule amount to the total halo mass, the 
only way to reveal the dark matter distribution is 
through dynamical analyses. For this purpose, the 
dynamically cold streams discussed in this paper 
are ideal. 
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Fig. 1. — The end-state of the simulation (at T — 12 Gyr) of the dwarf galaxy model Dl (originally containing 
10 5 particles) in the Galactic potential model HI with v c = 220 km s -1 , and q m = 0.9 is shown. The upper 
panel displays color-coded Heliocentric radial velocities, while the bottom panel shows the stellar distances, 
color-coded as a function of g* apparent magnitude. Though the distribution of the majority of the debris 
follows approximately a Great Circle on the sky, a significant amount of material is also precessed away from 
the Great Circle stream. The SDSS commissioning rvgi dataset is confined to ±1.25° of the equatorial plane 
in two stripes at 23. 5^ < RA < 3.5 h and 9.7 h < RA < 15. 5' 1 . The model predicts a high density of stream 
stars in the SDSS dataset at two regions at right ascension of ~ 2 hours and ~ 14 hours. This is where the 
A-colored star clumps are found in the SDSS dataset. 



a = 90 





Disruption time (Gyr before present epoch) The diagram shows a 2.5 de- 
gree wide slice through the (celestial) equatorial plane of the simulation displayed in Figure 1, which was 
evolved in a Galactic potential with an almost spherical halo (mass flattening q m = 0.9). The color of the 
particles marks the time since they were ejected beyond the tidal limit of the dwarf galaxy progenitor. The 
similarity of this distribution to the data of Yanny et al. (2000) (their Figures 3 and 11) is striking. In this 
model, the azimuthal ring-like features are produced by material that was disrupted off the Sagittarius dwarf 
galaxy at early times, and which has had a long time to precess in the Galactic potential. 
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a = 90 




Disruption time (Gyr before present epoch) Identical tO Figure 2, eXCept 

that the data is derived from a simulation in a q m = 0.5 flattened halo. Note the marked difference to 
Figure 2, especially in the right ascension range 270° < a < 0°. Further SDSS data in that region will 
greatly help the discrimination of competing models of the dynamical evolution of the Sagittarius dwarf 
galaxy, and of the structure of the outer Halo. The hole in the particle distribution near 210° is due to the 
structure of the tube orbits populated by the particles (see Figure 10 of paper 1). 
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